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Updated analysis of the dynamical relation between asteroid 
2003 EHl and comets C/1490 Yl and C/1385 Ul 
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ABSTRACT 

The asteroid 2003 EHl, proposed as the parent body of the Quadrantid meteor shower, is 
thought to be the remnant of a past cometary object, tentatively identified with the historical 
comets C/1490 Yl and C/1385 Ul. In the present work we use recovery astrometry to ex- 
tend the observed arc of 2003 EHl from 10 months to about 5 years, enough to exclude the 
proposed direct relationship of the asteroid with both of the comets. 

Key words: comets: individual (C/1490 Yl, C/1385 Ul) minor planets, asteroids: individual 
(2003 EHl) - astrometry - celestial mechanics - methods: numerical - methods: statistical 



1 INTRODUCTION 



' The Near Earth Object known as 2003 EHl was discovered in 
] March 2003 by B. S kiff in LONEOS survey data, and announced 
'. in MPEC 2003-E27 jMcCluskv etai]|2003l) . As first pointed out 
' by |jennisken3 ( l2004h . its orbit has a good similarity to that of the 

Quadrantid meteor stream, the last major annual meteor shower to 

be linked with a known comet or minor planet. 

Before this proposed linkage, many other asteroids and comets 
have been suggested as possible progenitors of the shower. Th e 

, most convincing of these objects, proposed bv lHasegawal i ll979l) . 

• is known as C/1490 Yl. Appearing as a bright naked-eye comet in 
the last days of 1490, it was extensively observed for the first two 
months of 149 1 by Chinese, K orean, and Japanese astronomers, 
as reported by IHo & Hoi ( Il962h . Various parabolic orbits for this 
object, based on the written descriptions of the comet path in the 
sky, have been proposed in the past, with the most recent given by 
iHasegawa (1979). 

There have also been various attempts to integrate the orbit o f 
the Quadrantid stream to 1491 (such as by lWilliams & Wu|[T993h . 
with the integrated orbit showing a remarkable similarity with 
the comet. They also proposed a possible identity of C/1490 Yl 
with another historical comet, C/1385 Ul, and used this possible 
linkage to estimate an eccentricity e 0.75 for C/1490 Yl. 

The comparison between various orbits is meaningful only if 
they are referred to the same epoch, because the action of various 
gravitational perturbations is often significant on a time span of 
a few centuries. This caution is especially important in this case, 
because the orbit of 2003 EHl has an aphelion around the orbit of 
Jupiter, therefore it often enters a region of space in which strong 
perturbing effects are acting. This comparison can only be done by 
integrating the orbit of 2003 EHl to 1491 (or 1385), the date of the 



perihelion of the comets, because the two cometary orbits are both 
too uncertain to be integrated forward meaningfully. 

Th e first attempt to perform such an integration to 1491 was 
done by I Williams et a l. (2004a, b), using observations of 2003 EHl 
spanning a time period of less than one year. They sampled the 
actual uncertainty region for 2003 EHl with 3500 synthetic clones, 
and integrated their orbits back to the 1491 epoch. The main 
result of their work is that the integrated position of the clones 
is extremely variable, but at least one of them had characteristics 
closely resembling those of comet C/1490 Yl, and its path in 
the sky also closely matches the one observed for the comet. 
However, most of the integrated orbits predicted the object to be 
around aphelion, and not close to the perihelion at which comet 
C/1490 Yl was observed. They didn't attempt any integration to 
1385, to compare the orbit with C/1385 Ul, but they used the 
eccentricity inferred from the identification in their comparisons. 

In April 2008 we recovered 2003 EHl using the University of 
Hawaii 2.2 m telescope on Mauna Kea, Hawaii, equipped with a 
Tektronix 2048x2048 CCD camera. We observed it for two nights 
and extracted 7 accurate astrometric positions from our set of 
images. These data extend the observed arc by a factor of more 
than 5 and should allow a comparable reduction in the expected 
positional uncertainty for the orbital extrapolation to 1491. The 
result of this updated orbital comparison, referred to both the 1491 
and 1385 comets, is presented in this work. 



2 ANALYSIS 
2.1 Procedure 

In our analysis w e used a method similar to the one applied 
by IWilliams et al. I ( l2004ij) . We first generated 1200 stochastic 
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clones of 2003 EHl and integrated their orbits to January 1491. 
We repeated the same approach, creating 600 other clones, and 
integrated them to 1385. The main advantage of this approach is 
that it doesn't require any linear approximation, and it is therefore 
suitable considering the large uncertainties associated with this 
situation. 

The orbit generation procedure was done using the Monte 
Carlo feature of the software "FindOrb" written by B. Gra}!]; 
the procedure creates the stochastic orbits by adding Gaussian 
noise of given amplitude to each astrometric observation, and 
then used this new set of simulated positions to compute a best-fit 
orbit. These orbits were then integrated back to the desired epoch, 
1491 January 22 in our case (the same mid-arc time used in 
IWilliamsetal]|2004ah and 1385 November 01 (the perihelion 
date of C/1385 Ul). We used a Gaussian noise model with 0.8" 
as the standard deviatior0. The integration takes into account the 
gravitational influence of all the major planets, the Moon and Pluto 
according to the long-term planetary ephemeris JPL DE-406, and 
of (1) Ceres, (2) Pallas and (4) Vesta. 



A similar work applied to the 1385 apparition gives a compa- 
rable disagreement: all the clones are spread on a line of variation 
~ 100° in length and centered at ~ 100° from the expected position 
of C/1385 Ul. The closest clone is still at about 65° from the 
comet position. 

The present orbital elements of 2003 EHl (and their cor- 
relation matrix) are shown in Table [T] while the integrated ones 
are compared with those of the two historical comets in Table [2] 
According to the nominal orbit, during the integration interval the 
object had close encounters with Venus (within 0.04 ua in 1517), 
Earth (within 0.10 ua in 1818) and Jupiter (within 0.25 ua in 1664), 
all capable of significant orbital perturbations. 

Although our new observations seem to exclude the identifi- 
cation of 2003 EHl with C/1490 Yl and C/1385 Ul, we cannot 
exclude the possibility that 2003 EHl and the comets are fragments 
of some parent body that split long ago. 



2.2 Results 

The evaluation of the resul ts has been done in a way consistent 
with lWilUams et al.l j2004ah ; we used the stochastically generated 
orbits to compute the coordinates of each clone in the plane of the 
sky and compared them with the position described in the ancient 
records. 

As e xpected, our cl o nes lay on the same line of variation 
found by Iwilliams et al] ( l2004ij) . but the spread of the points 
along the line is ~ 5 times smaller, being ~ 20° and mainly in 
declination. On the reference date of 1491 January 22 the region 
is centered around 6 - -58°, and the northernmost clone is at 
6 - -46°, at about 85° fr om the reported positio n of the comet. 
As already pointed out by IWilliams et alj ( l2004 al) these positions 
correspond to objects close to the aphelion of the cometary orbit. 
The path in the sky is completely different from the cometary one, 
and the velocity of the clones is lower, as expected knowing that 
most of them are located around aphelion. None of the clones is 
consequently compatible with the information given by the ancient 
records. 

The computation described above neglects the possibility of 
non-gravitational effects acting on the object. Under conservative 
assumption!^ we estimate that the northernmost clone will still 
need a non-gravitational acceleration A > 2 x 10"** ua d"^ to 
have a perihelion date compatible with C/1490 Yl. This value is 
much larger than a typical cometary non-gravitational parameter, 
and comparable only with the ones measured for the most active 
comets known so far. 



' Downloadable at http : //www . pro j ectpluto . com/£ind_orb . htm] 
^ This is probably a strong overestimation of the true astrometric error of 
our Mauna Kea observations, which are internally consistent to the ~ 0.2" 
level. 

^ We assumed that all the non-gravitational acceleration acted on the mean 
motion, that the average gravitational acceleration was equal to its value 
at aphelion and that the object was continuously active from 1491 to the 
present. 



3 CONCLUSIONS 

The results obtained in this work seem to exclude the proposed 
identification of comet C/1490 Yl and C/1385 Ul as historical 
cometary apparitions of the asteroid 2003 EHl. None of the in- 
tegrated orbits matches the position and the motion of the 1491 and 
1385 objects. 

The identity between the two historical comets obviously 
cannot be excluded by this analysis. Similarly, the possible identi- 
fication of 2003 EHl as the progenitor of the Quadrantid meteor 
shower is still plausible, although it now lacks of any confirmed 
cometary activity from the parent body. 
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Table 1. Orbital elements and correlation matrix for 2003 EHl at the ciurent standard epoch (JD 2454600.5) and mean ecUptic and equinox (J2000.0). 
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Table 2. Integrated orbital elements of the asteroid compared with the two historical comets (mean ecliptic and equinox of J2000.0). The error bars are the 
99.7 % limits for the distribution of each element. Most of the distributions are strongly non-Gaussian, and some of them are bimodal, presumably due to close 
planetary encoimters. 



